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Abstract 

We consider the problem of an ideal polymer confined in a droplet. When the 

droplet radius is smaller than the (unconfined) polymer radius of gyration, the 

polymer entropy will depend on the droplet shape. We compute the resulting 

surface free energy. Using parameters appropriate for polymers confined in 

microemulsions, we find that the polymer and bending surface energies are 

comparable for the lowest modes. Finally, we argue that chain self-avoidance 

will decrease the strength of the polymer contribution to the surface energy. 
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The Hamiltonian governing surface shapes usually follows from considerations of two- 
dimensional elasticity and geometry Thus, fluid surfaces such as bilayers, emulsions and 
microemulsions are well-described by surface tension and curvature elasticity, while solid 
membranes also depend on in-plane shear and compression moduli jl]. In general, however, 
one must also take into account the response of the surrounding three-dimensional medium. 
(For convenience we speak of only one medium although in general the surface may separate 
two different materials.) The medium is characterized by a length, £, which describes the 
distance over which disturbances are screened in the bulk. If we consider only surface 
perturbations with wavelengths larger than £, the bulk will contribute to the surface elastic 
constants but the form of the Hamiltonian will be as described above. For perturbations 
with wavelengths below £, on the other hand, the bulk response gives rise to new terms in 
the surface Hamiltonian. 

As an illustration, consider a surface bounding an elastic solid with modulus B. In this 
case, £ is infinite. The contribution to the surface energy from the bulk elasticity can be 
estimated by dimensional analysis for surface perturbations with wavelengths much less than 
the size of the system. If we allow a surface tension, /i, and bending rigidity, k, then the 
energy per area of a surface displacement u(q) exp(igr) must take the form: 

Energy/Area ~ (Bq + fiq 2 + Kq 4 )\u(q)\ 2 . (1) 

Thus, for sufficiently small q, the bulk elasticity dominates. 

In this paper we consider an approximately spherical surface (droplet) enclosing a single 
polymer. We are interested in the limit in which the radius of the droplet is smaller than the 
radius of gyration of the unconfined polymer. We were motivated by recent work in which 
individual polymers have been confined within microemulsion droplets 0, 0. In order to 
make the problem tractable, we consider ideal polymers, i.e. we neglect oidance. We will 
return to this point at the end. 

The droplet is assumed to be incompressible, with volume V = 4irR 3 /3. Points on the 
droplet surface, S, are parametrized in terms of the radial displacement from the sphere 
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with radius R: 



Incompressibility implies: 



Rs(6, ( j ) )^(R + u(9, ( p))r. (2) 



dn^R s (e,<f>) 3 = ^. (3) 



After expanding u in normalized spherical harmonics 



oo I 

«(M)=E E u lm Y lm (9,(f>), (4) 

Z=0 m=-Z 

we enforce the constraint Eq. (|3|) by solving for uoo'- 

Uoo = -^^Y,\ u im\ 2 + 0{ufj. (5) 

The surface free energy of the droplet has the form JF = jF b + jF p , where T v is the free 
energy of the confined polymer and jF b is the bending energy of the droplet surface (c.f. ||): 

Th= \ E %^ + 2)(i - mil + l)-4w + 2w*). (6) 

i>0,m 

K is the bending rigidity and w is the ratio of R to the spontanteous curvature radius. For 
simplicity, we have ignored surface tension, although it may easily be included. T v is the 
free energy of the confined polymer. 

For an ideal polymer of N segments, T v may be expressed in terms of the statistical 
weight Gjsf(ro,ri), which is the number of walks that start at r and end at r\ after N steps: 



T p = -k B T log [j v d 3 r J^n G?iv(ro,ri)J . (7) 

Gjv( r O; r i) satisfies the diffusion equation (c.f. [|]) 

— G N (r , n) = -V^G^ro, n), (8) 

where b is the effective bond length. The radius of gyration of the polymer, R g , is related 
to iV by Rg = Nb 2 /6. For either r or r x on the boundary of the droplet, S, Gjv(ro,ri) 
vanishes. In terms of orthonormal eigenf unctions ipi of the Laplacian, satisfying 



ViJ)i = -kfy i ^(r)| r6S = 0, (9) 

we have the expansion 

G N {r , ri ) =^V*(r )Vi(r 1 )exp(-^^). (10) 

i 

Since S is a perturbation of a sphere of radius R, Eq. (||), the eigenvalues k 2 are proportional 
to 1/R 2 . In the limit i? 3> Crjv(ro, i"i) falls off rapidly for |r — r x | R g , and Eq. (^) 
reduces to the translational free energy of the polymer, as expected: T v ~ —k B T log(V). 

For R g > R, the sum in Eq. fllCf ) can be approximated by keeping only the lowest 
eigenvalue k 2 (ground-state dominance ||). Within this approximation, the free energy is 
given by 

T v = k B TR 2 g kl - 2* fl Tlog(c b ) + O ( e -^( fc ?- fc o)) 
c = f d 3 r^o(r). (11) 

The first term on the right-hand side of Eq. ([□]) is the leading term (in R/R g ). We have 
also included the first correction (the second term in Eq. (|TT|)) because it may be calculated 
with little extra work (see appendix). Note that the remaining higher order corrections 
are exponentially suppressed. Thus, for a confined chain, the shape dependence of T v is 
determined by the shift in k\ and ipo for perturbations of the form Eq. (|2|). 

As a simple application of the above formula, we first compute the free energy of the 
polymer as a function of the droplet radius. We therefore take Rs = R and the eigenvalues 
in Eq. @ reduce to (mr/R) 2 , n = 1, . . .. Thus, for R < R g , the free energy is given by 

T v « k B Tn 2 ^ 2 . (12) 

It is tempting to view the interior of the droplet as an elastic solid with bulk modulus B. 
To estimate B, we relax the requirement of incompressibility and compute the change in free 
energy under a uniform change in the radius of the droplet. We thus take R —>■ R — u and 
expand Eq. (|T^) in the strain u/R; the term linear in u/R simply reflects the fact that we are 
perturbing about a compressed chain. Dividing by the droplet volume gives a bulk modulus 
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of order B ~ ksT R^/R 5 . We can now use the result from the introduction to estimate T v 
for surface disturbances with wavelength much smaller than R. This corresponds to a mode 
ui m with wavevector q ~ l/R, I ^> 1. From the argument in the introduction we then find 



T ^ kBTl %%- (13) 



As we shall see below, T v does indeed take this form in the limit of large /. However, it 
is in the low q modes where the contribution from the polymer is expected to be important 
(compared with the contribution from the bending energy). This corresponds to small / (i.e. 
wavelengths comparable to R) and requires a more careful analysis, which we now present. 

As described above, we must find the ground state of Eq. (|9]). The general problem of 
the change in the spectrum of the Laplacian under a boundary perturbation was considered 
some time ago and appears in a number of contexts, such as acoustics in or around irregular 
regions, cavity resonators, and the Wigner-Seitz approximation for electrons in a crystal 
(see ]7[ and references therein). The early work on this problem focused on the solution to 
first order in the perturbation. Unfortunately, we require the result to second order (i.e. 
0(uf m )), which is substantially more difficult to determine. In ]7| a general formula is given 
for the change in the eigenvalues to second order, however the final formula appears to be 
incorrect. Therefore, in the appendix we calculate the ground-state eigenvalue and integral 
of the ground-state wavefunction for a perturbation of the form Eq. (||). 

Combining the results of the appendix with the constraint Eq. ([5]) gives: 

^[i+E(^-f)^ + o(<4)] (i4) 

t>0,m 



'2 



7T 



c = 2^R^ 1 + ~ £ {2(1 + 2)/, - TT/f - 3 ^)^! + o(uf m : 



l>0,m 

where fi = ji+iijt) / and ji is the spherical Bessel function of order /. Note that the 
Ui m -dependent terms in the volume preserving perturbation Eqs. (0,§) correspond to in- 
finitesimal translations, which should not change the spectrum of the Laplacian. One may 
readily check that indeed the order u\ m terms in Eqs. ( JT^ - P^D vanish. 
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Substituting the above results into Eq. ([□]), and dropping an additive constant that is 
independent of U[ m , we find for the polymer contribution to the droplet free energy: 



|2 



^ = ^£%r I + ... (is) 



r 



2 ,6? R 2 

(>0,m 

Rl „ „ ,x 7T „, Z + 2 , 3/ + 6 



+ 2 - 7T/,) + t# - + 



In the limit of large Z, T; has the behavior 

r,^M§ + §). (is) 

We thus recover the form predicted from simple elasticity Eq. (P~3"D. 

As discussed above, the I = mode is constrained by incompressibility and the I = 1 
modes correspond to translations (which implies I\ = 0). Therefore, the leading contribution 
in the small-/ or long- wavelength limit comes from 1 = 2: 

We thus find that for R g > R, the surface free energy picks up a contribution from the 
polymer which does not appear to have a simple interpretation in terms of surface elasticity 
and geometry. Furthermore, for the lowest modes, this contribution to the surface energy 
can be quite important. For the experiments in 0, R g /R ~ 4, k w 2&bT, and w « 2. With 
these values, the ratio of the (I = 2) polymer contribution to the free energy to the bending 
contribution is ^{1=2/^1=2 ~ 2.4; the contributions from bending energy and the polymer 
are comparable. 

In the calculations above we have neglected self- avoidance. However, in the limit of very 
high compressions self-avoidance will certainly be important. R is bounded from below by 
R m , the radius for which the polymer volume fraction is equal to one. Near this limit, we 
expect the interior of the droplet to be similar to a melt. The screening length £ is then 
microscopic and the polymer free energy is insensitive to the shape of the droplet: w 0. 
Therefore, for the polymer to make a significant contribution to the surface free energy, the 
droplet radius must be in the range R m < R < R g . However, even for droplets within this 



range, self avoidance will be important for the lowest modes. A simple scaling estimate of £ 
§ gives 

f- fl U) • (i8) 

where R g ~ N u . For a chain in a good solvent, z/ « 3/5, and in a theta solvent, v = 1/2 
f|. Thus, even for a modest compression of the polymer (i? < the screening length is 
comparable to R (£ ^ i?) and the lowest modes of the droplet are just at the wavelengths 
where screening becomes important. While this results in a suppression of the free energy 
of the lowest modes, we cannot gauge the size of the effect. In particular, we do not know 
whether this suppression is strong enough to significantly reduce the polymer contribution 
to the surface energy relative to the bending energy. 

Even for the higher modes, for which screening is unimportant (i.e. I > R/£), there will 
be corrections due to self- avoidance. We can again use a scaling argument to estimate the 
free energy of a polymer in a spherical droplet (with R < R g ): 

fR \ l/v 

?P*k B T(f±) . (19) 

Ignoring numerical factors, the above estimate is smaller than the ideal polymer result, Eq. 
(|12|), by a factor of (R g / R) 2 ^ 1 ^ . Following the arguments leading to Eq. fljjD , we then 
expect the large-/ behavior of T/ to be smaller by a factor of (Rg/R) 1 ^ 3 for a polymer in a 
good solvent; there is no correction for a polymer in a theta solvent. 

To summarize, we have determined the contribution to the surface free energy arising 
from confining an ideal polymer within a droplet. Not surprisingly, this contribution does 
not have a simple interpretation in terms of two-dimensional elasticity and geometry. Using 
parameters appropriate for polymers confined in microemulsions we find that the bending 
and polymer terms in the surface energy are comparable for the lowest modes. Scaling 
arguments suggest that chain self-avoidance will tend to decrease the magnitude of the 
polymer contribution. 
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APPENDIX: BOUNDARY PERTURBATIONS 



In this appendix we find the ground-state eigenvalue and integral of the ground-state 
eigenfunction to second order in the perturbation of the shape of the boundary, Eq. (0). In 
a formula for the eigenvalues of the Laplacian to second order in a boundary perturbation 
is given for the case of Dirichlet boundary conditions (Eq. (9.2.71) of |7|]). We have not 
been able to reproduce the derivation of this formula. Furthermore, when one applies the 
formula of [[JJ to the simple case of a spherical boundary with a perturbation of the radius, 
R — > R — 5R, it does not give the correct result. However, by a similar procedure to that 
described in [0, we can derive the correct formula for general perturbations of a sphere. 

We follow the notation in [|7| and denote the perturbed volume by V and the boundary 
by S. S is a perturbation of So, which is a sphere of radius R. For the purposes of the 
derivation, however, we must formally consider perturbations that are entirely contained 
within the unperturbed surface. We therefore consider a second sphere S' with a slightly 
larger radius than So, such that S is entirely contained in S' (Fig. [I]); the volume contained 
in S' is V. We will compute the spectrum for the Laplacian on V in terms of the spectrum 
on V . The surface S is thus given by: 

T S = (B! + v(6,4>))i u(M)<0, (Al) 

where in terms of the perturbation u(9, 0) ( Eq. (0)) 

v(e,(f))=u{e,(j))-a R' = R + a 

a = max it(0, 0). (A2) 

9,(f> 

We wish to determine the ground-state eigenfunction and eigenvalue on V: 

(V 2 + fc 2 )</> = ^(r)| re5 = 0. (A3) 
in terms of the orthonormal eigenfunctions and eigenvalues on V: 

(V 2 + k 2 nl )<p n i m = 0nim(r)\res' = 0. (A4) 



Since S' is a sphere of radius R', Eq. ([A4]) is easily solved: 



*m = ^ (A5) 

where Z = 0, 1, . . . , m— —I, . . . I, n = 1, 2, . . . , ji(x) is the spherical Bessel function of order 
I and x n i is the nth root of ji(x) ||. To leading order we have ip(r) = 0ioo( r ) + 0(v) and 
k 2 = kf Q + 0(v ). A simple application of Green's theorem gives: 

(k 2 - kl ) J y d 3 r loo (r)^(r) = - jf d 2 r loo (r) (|^^ , (A6) 

where we have used the fact that ?/>(r) vanishes on S; d/dn denotes the outward normal 
derivative at the boundary. As in [0 , we solve this integral equation iteratively in powers of 
v; the right-hand side is 0(v), since 0ioo( r ) vanishes on 5". 

At the lowest order, we substitute 0ioo(j) for ip{r) in Eq. (|A6[) , express the integrals over 
V and S in terms of the nearby V and S', and expand in v to give 

k 2 - A;2 Q = — l—k l0 hi v lm + 0(v 2 ) , (A7) 
2R W n hn 

where V\ m are defined by expanding v(9, (f>) in spherical harmonics: 

v{e,<f>) = y £vi m Yi m {e,<i>). (as) 

lm 

By a further application of Green's theorem, and again using Eq. ( |A6| ) , the ground-state 
eigenfunction may be written as (see Rj) 

ip(r) = 0ioo(r) + J'V E'^ d^^ 0n;m(r) + 0(w 2 ). (A9) 



The prime on the sum in Eq. (|A9|) indicates that the state (n = 1, / = m = 0) is excluded. 
In order to compute k 2 to second order in v, we would like to substitute the above expression 
for ip(r) into Eq. (|A6|) . Unfortunately, the sum over n in Eq. ( |A9|) converges quite slowly 
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(k n i ~ n for large n) and we cannot move the derivative in Eq. ( |A6|) past this sum. In order 
to improve the convergence, we sum the leading large- n behavior in Eq. (|M|) usingQ 



4>nlm{r) 



VlmXl 



I'm 



nlrn lhn ^ 

which follows from the identity (for r < 1 



£ v2/?' V/?' 



(A10) 



n^X %nljl+x{Xnl) 



E 



r 
2 



(All) 



After summing the large-n behavior, Eq. (|A9|) gives: 

k 



i>(r) = (j>xoo(r) + 



40 



k 2 



,[E^ m i, ( ,2 " t, 2 \<t>nlm{r) 



(A12) 



nlm 



+ E 



vi m Y lm (6, 4>) ( r 



0ioo (r) 



k 
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+ 0(v 2 ). 



lm V2R 

This expression converges sufficiently rapidly that it may be substituted into Eq. 
We evaluate the resulting sums over x n \ by contour integration || and then use Eq. (|A2|) to 
express the perturbation in terms of u. The final result for the ground state eigenvalue is: 



k 2 - — 
k ~ R 2 



3 u 



\Ulr 



R 2 



-(l + l-irfl. 



0(u 3 



(A13) 



V^R ^R 2 

where ft = 

By a similar iteration, again using Eq. ( |A10| ) to improve convergence, we find for the 
integral of the normalized ground-state eigenf unction: 



fv d 3 r^(r) 
f v d 3 r V>( r ) 2 



U-R?' 2 



7T 



1 + 



3 Mqo 3 u 2 



+ 



00 



4^ R 32tt R 2 



(A14) 



+ \ £ ( 2 (i + 2)/ I -7r/f-3- ( ' +l: 

l>0,m 



71 



R 2 



1 We note that the identity Eq. ( |A10| ) is a counterexample to a claimed general result used in JjJ 
to improve the convergence of sums arising in boundary perturbations (see p. 1044 of (j). 
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FIGURES 

FIG. 1. The surface of the droplet, S, is a perturbation of a sphere So- S may be viewed as a 
perturbation of a slightly larger sphere, S', which contains S. 
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